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Abstract 

We study the Abrikosov-Nielsen-Olesen multi-vortices. Using a numerical code 
we are able to solve the vortex equations with winding number up to n = 25, 000. We 
can thus check the wall vortex conjecture previously made in pq. The numerical 
results show a remarkable agreement with the theoretical predictions. 
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1 Introduction 

The idea of the wall vortex emerged in P] when we faced the following situation. 
Consider a theory that has two degenerate vacua, one in the Higgs phase and another 
in the Coulomb phase. Then take a domain wall that interpolates between the two 
vacua and place a monopole on the Coulomb side. What happens if we continuously 
move the monopole towards the Higgs phase? It is known that a monopole in the 
Higgs phase cannot exist by itself and must be confined. The only reasonable thing 
that can happen is that the monopole wall system is continuously transformed into a 
monopole-vortex-wall. This implies that the confining vortex in the Higgs phase can 
be continuously transformed into a wall or, in another way, is made of the same stuff 
as the wall. Thus we call it a wall vortex. Now it is easy to understand what are the 
forces that keep the vortex soliton together. The wall tension that tends to squeeze 
the vortex, is nothing but the Derrick j3] collapse force coming from the scalar part 
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of the action d(j)d(j) + V{(j)). What prevents the soliton from collapsing is a pressure 
coming from the magnetic flux inside the vortex. The magnetic energy is but 

we should keep in mind that when varying the radius of the cylinder, what remains 
constant is not the magnetic held but the hux. In terms of the hux, the magnetic 
energy is ^b /and thus a pressure term. 

The wall vortex argument, up to this point, was only a qualitative way to under¬ 
stand the continuous transition from a wall to a vortex and also the forces that bind 
together the soliton. In order to apply this idea also quantitatively we need to hnd a 
regime in which the radius of the vortex Ry becomes much greater than the thickness 
of the domain wall Ayr. We hnd that this is exactly the large n limit where n is the 
number of quanta of magnetic huxes carried by the vortex. It is easy to understand 
why this happens. The domain wall does not know anything about the magnetic hux 
and thus its thickness is independent of n. On the contrary the radius of the vortex 
depends on n, in particular it grows as n grows. This means that the ratio Ay/ / Ry 
can be made arbitrarily small by increasing the parameter n. 

In j2] we applied the wall vortex argument to the general Abelian-Higgs model. 
In this case there is only one true vacuum, the Higgs vacuum. The Coulomb phase 
is not a true vacuum of the theory, but due to the symmetry of the problem, it is 
always a stationary point of the potential. Is it possible, that also in this case, in the 
large n limit, the vortex becomes a wall vortex? This was the conjecture made in j2]. 
The aim of this paper is to prove this conjecture by means of numerical computation. 

We organize the paper in the following way. In Section |21 we give a review of the 
wall vortex idea and we state the theoretical arguments that support the conjecture. 
Then in Section IHl we present the numerical proof of the conjecture. In Section Id.21 
we analyze the original problem when the Coulomb and Higgs phases are degenerate. 
Finally in Section 0] we conclude with a discussion about the relation between different 
potentials in the large n limit. 

2 Theoretical Section 

In this section we provide the theoretical analysis that leads to the wall vortex 
conjecture. A lot of independent arguments support the wall vortex limit. Part of 
these arguments are new and part are just a review of |I1I2|. 
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The theory under consideration is the Abelian-Higgs model, which is the relativis¬ 
tic version of the Ginzburg-Landau theory of superconductivity ji]. It is a U{1) gauge 
theory coupled to a charged scalar field q 

C. = -4^^’,--^'“' - 1 ( 3 ^ - - r(l9l). (2.1) 

The potential G(|g|) is such that the field q acquires a vev go and gives mass to 
the U{1) photon. In this phase the theory admits vortex solutions called Abrikosov- 
Nielsen-Olesen (ANO) flux tubes |S1 El- The ANO vortex is a soliton extended in 
1-1-1 dimensions. The usual way is to choose cylindrical coordinates {z, r, 9) with the 
vortex oriented in the i direction. The fields can then be put in the following form 

q = , ( 2 . 2 ) 

71 

Ae = -A(r), 

r 

where g(r) and A{r) are profile functions that must be determined using the equation 
of motion. The field q at r ^ 00 has been chosen to lie in the vacuum manifolds 
|g| = go. This is necessary to have a finite energy configuration. The solution ()2.2|1 
is obtained by choosing the element n of the homotopy group [n] G 7ri(S^). To have 
finite energy it is also necessary to turn on a gauge field in order for the covariant 
derivative Dq to vanish. This creates a magnetic flux that is exactly proportional to 
the topological number n. 

2.1 The Wall Vortex Limit 

Now we describe what is the wall vortex and then we claim that multi-vortices are 
wall vortices in the large n limit. 

The main idea is that the Coulomb phase g = 0, even if it is not the true vacuum 
of the theory, is always a stationary point of the potential. This in fact is a trivial 
observation but, as we will see, it is full of non-trivial consequences. In general there 
are two cases described in Figure [H The Coulomb phase is a false vacuum with energy 
density V (0) = £0 and can be metastable or instable. We call Vq the highest peak of 
the potential between the Coulomb and the Higgs phase. 

Now comes the first assumption. Consider the domain wall interpolating between 
the true and the false vacuum as an independent object with tension and thickness. 
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Figure 1: Two possibilities for the Higgs potential. In the first case the Coulomb phase is instable 
while in the second case it is metastable. 

respectively: 

Tw ~ 5 ~ ■ (2.3) 

This domain wall does not really exist as an independent object. The energy density 
Eq gives a negative pressure that tends to eliminate the Coulomb half space leaving 
only the Higgs vacuum. But, as we will see, it can exist in the vortex solution where 
the negative pressure is counter balanced by the pressure of the magnetic held. The 
key point is that the domain wall maintain its “identity” even if we change the quanta 
of magnetic hux n. Thus in the large n limit the radius of the vortex Ry becomes 
large while the thickness of the wall remains hxed. This is the key assumption 
that needs to be proven. 

For the moment let’s assume that the wall has its own identity and then go on to 
the consequences. The wall vortex is a hux tube obtained in the following way. We 
compactify the domain wall on a cylinder of radius i?, keeping the Coulomb phase 
inside and the Higgs phase outside. Then we turn on a magnetic hux inside the 
cylinder and we write the tension of the tube as function of the radius : 

27m^ 

T{R) = ^2j^2 T\y2'iiR + EqttR^ . (2-4) 

There are three forces that enter the game. Two of them tend to squeeze the tube 
and are due to the tension of the wall and to the energy density ^o- This is nothing 
but the Derrick collapse coming from the scalar part of the Lagrangian. The last 
force is a positive pressure that comes from the magnetic held inside the tube. 
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There are two regimes in which (E3D easily can be solved: the SLAG bag and the 
MIT bagd 


SLAG bag: This regime is when n satisfies the two conditions 


qo^evo 

■C n -C 




^0 


3/2 


(2.5) 


The first condition is that the radius Rv is much bigger than the thickness of 
the wall Aw- The second condition is that the surface term Tw^ttR dominates 
over the volume term EqTtR^. In this limit 
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SLAG 


= 35 ^- 


TT 


rp \ 2/3 
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-Rslac — ^ 




n 


2/3 


( 2 . 6 ) 


MIT bag: This regime is when n satisfies the condition 


qo^evQ 


^0 


3/2 


<C n . 


(2.7) 


In this limit the volume term in (E3D dominates over the surface term and 


^MIT — - n 


R 


MIT 


= ^- 


e 

Note that the tension is proportional to n, as in the BPS case. 


n . 


( 2 . 8 ) 


2 

Let’s write for clarity the complete story of multi-vortices. At the value n S> 

where the radius is much bigger than the thickness of the wall, the multi-vortex 

becomes a wall vortex. If the parameter uq is much bigger than £o; the wall vortex 

2 

can be subdivided into two different regimes. In the first one n <C and the 

surface term dominates over the volume term and we call it the SLAG bag regime. In 
this regime the tension scales like When n ~ ; fhe volume term starts to be 

comparable to the surface term and thus a second order phase transition between the 
SLAG bag and the MIT bag takes place. Note that the MIT bag regime, whenever 
£o 7 ^ 0, is always present and always dominates in the large n limit. 

Now lets write the conjecture for clarity: 

Consider the Abelian Higgs model (Unp with a general potential that has a true 
vacuum at \q\ = q^ ^ 0 and a Coulomb phase with energy density V{0) = Eq ^ 0. 

^We have taken the names from the bags models of hadrons: the MIT bag [H] and the SLAG bag 
cni In the Friedberg-Lee model m they arise as different limits. 
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Call Ty{n) the tension of the vortex with n units of magnetic flux. The conjecture is 
that 


lim Tv{n) = TuiTin) . (2.9) 


2.2 The Differential Equations 


The differential equations for the prohle functions of the vortex ( 1221 ) are 

^ + 2 (1 - = n 

r dr ^ ^ 2 Sq ’ 


( 2 . 10 ) 


d^A 1 dA 


+ 26^(1 - A)q^ = 0 , 


r dr 

where n is the winding number. We are looking for some limit of the parameters so 
that the vortex really looks like a wrapped wall. In this limit the prohle functions 
should be: 


qir) = qo9H{r-Rv), ( 2 . 11 ) 

A{r) = r'^/Ry 0 < r < Ry , 

A{r) = 1 r > Ry , 


where 6h is the Heaviside step function. Note that the magnetic held is oc Aflr and 
is constant inside the vortex and zero outside. 


First of all we manipulate the diherential equations (ITTUl) to simplify them. The 
potential can be written as a dimensionless function 


!/(«) = £„V 



( 2 . 12 ) 


where £0 is the value of the potential at g = 0 and go is the vev at which the potential 
vanishes. The rescaled potential fullhlls V(0) = 1 and V(l) = 0. We also rescale the 
scalar held g = go X- In the case of a quartic potential the only one with the required 
properties is V(x) = (y^ — 1)^- After these rescalings the equations (I2.1()j) for the 
prohles become: 


d‘^X 

dr"^ 

d^A 

dr‘^ 


r dr 


1 dA 
r dr 


+ b{l-A)x^ 


dV 

X - = 0 , 

W 


= 0 . 


(2.13) 

(2.14) 
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There are three parameters that enter the game: 


n , 


1 ^0 





h = 2e^go^ • 


(2.15) 


Now comes the first non-trivial hint for the wall vortex conjecture. If the wall limit 
exists, then formula (EH) can be trusted in this limit. But the radius of the vortex Ry 
comes from equations (ITTUD and must depend only on the three relevant parameters 
n, a, b. In general a function of n, e, Eq, qo cannot be expressed as a function of n, a, b, 
but for (IZHD this is possible: 

= (2.16) 
Vab 

If we would not have found such an expression, we would have concluded that the 
wall vortex limit does not exist. This result encourages us to go on. 

In this paragraph we look for a limit in which the radius of the vortex remains 
constant and the solution approaches the wall vortex (ITTTl) . To do so we need to 
rescale also the parameters a and b such that the radius remains hxed while approach¬ 
ing the large n limit. The y prohle must become a step function: it is zero inside 
Rv, it goes from zero to one in a distance Aw, and then remains constant at one. 
Thus x" in equation develops a 5'(r — Ry) singularity or, in terms of Aw, 

x" ~ 1/Aw^- To counter balance this divergence in (I2.1dj) we must have that a goes 
to inhnity like 1 /Avk^- Now consider the second equation ()2.14|) where A"{Ry) has 
a {2/Ry)5{r — Ry) singularity, or in terms of Aw, A"{Ry) ~ l/{RyAw)- Since 
(1 — y4)x is of order Aw/Rv around Ry^ we must also send b to inhnity like 1/Aw^- 
Now we can reformulate the conjecture of the wall limit in another equivalent way. 

Consider the succession of parameters a„ = nai and bn = nbi and name the 
solution of \2.RA) and Ii 2 . 14 }) with the vortex boundary conditions, Xn,a„,6„(’") and 
An,an,bA^) ■ In the limit n ^ oo 


lim Xn,an,bAr) 

^^oo 

lim An,a„,bAn) 


nir - R 


-vj 


r 


jRy‘^ 0 < r < Ry 


(2.17) 


This limit has been chosen such that the radius of the vortex remains constant 
Ry = ^2/(ai6i) and also the ratio On/bn remains constant. The information about 
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the ratio a/b disappears in the wall vortex limit. It is only related to the shape of 
the limiting functions Xn,a„,bni^) ^^^1 Probably a stronger version of the 

conjecture is true: the ratio a„/6„ is kept limited from above and from below during 
the limit, so that it does not go neither to infinity nor to zero. 

The ratio a/b has also another important meaning. In fact, in the case of quartic 
potential, it is essentially the parameter (3 that measures the ratio between the Higgs 
and the photon mass: 

P = !^ = Sj. ( 2 . 18 ) 

b 

When /3 < 1 the Higgs attraction dominates and the vortices are of Type I. When 
(3 > 1 the photon repulsion dominates and the vortices are of Type H. When (3 = 1 
the vortices are BPS, this means that the tension is a linear function of n and the 
distance between different vortices is the vortex moduli space. 

2.3 A Non Trivial Check 

Now we make a non-trivial check of the result (El) using the famous example solved 
by Bogomol’nyi |I2|. When the potential is 

nkl) = y(l0^-O^ (2.19) 

the tension is 

Peps = , (2.20) 

for all n. Solving the model with our trick, the result must coincide with eq. 

For the BPS potential (EH, the energy density of the instable Coulomb vacuum 
is Eq = ^d‘2. Using ()2.8|1 . we hnd exactly (EH- This could hardly be just a 
coincidence. 

It is well known that n BPS vortices have a moduli space of real dimension 2n 
where the coordinates can be interpreted as the position of every constituent 1-vortex 
na. In the BPS bag we expect to recover this moduli space and we expect it to have 
infinite dimension. In fact this space consist of the closed surfaces with fixed area. 
More details will be discussed in ^3] in relation with the multi-monopole moduli 


space. 





3 Numerical Analysis 


This is the central part of the paper. With the help of numerical computations 
we are going to prove the conjecture made in the previous section. The strategy is 
to prove the conjecture as it has been reformulated in (EUD, that is, we rescale the 
parameters of the theory such that the radius of the vortex remains constant and the 
thickness of the wall goes to zero. 

The code^ used to solve the system of differential equations (j2.18112.14jl uses a 
numerical method called finite difference method. The boundary conditions are at 
the singular points of the system 0 and oo and thus we use the limits: 


lim {2A{r) — rW) = 0 , lim (nx(r) — rx) = 0 , 


lim 

r—j'CxD 




1 , 


lim 

r^oo 



1 . 


(3.1) 


Other methods often used are e.g. the shooting method, but in case of a step function 
insanely high accuracy is needed. Still for the finite difference method, the step 
function poses severe problems and the technique to deal with the problem is to use a 
combination of a continuation method and feeding the algorithm with an approximate 
solution. Also the rescaling of the equations such that the radius of the vortex is 
constant is an important factor for solving the equations. The finite difference method 
transforms the problem into a matrix equation 


Fij = 0 , i = 1,..., meshsize , j = 1, 2 , (3.2) 

which, when using a continuation method, looks like 

G,j{t) = 0, te[0,l], (3.3) 

where the problem is easy to solve when t = 0 and hard to solve for t = 1. Then 
the solver will try to increase the hardness in sufficiently small steps until the final 
problem t = 1 is solved. Whenever a solution has been found it is fed to the solver 
as an approximate solution which is the solution to the problem for f = 0 and for 
sufficiently small steps it is then possible to achieve higher winding number n. The 
approximate solutions tell the solver where to look for the solution in phase space. 

^The code is written in the mathematical language Maple. 
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3.1 Quartic Potential and the MIT bag regime 


The simplest case is the quartic potential. After the rescaling (ITT^ there is only one 
quartic potential to be considered: 

V(x) = - 1)2 (3.4) 


The following hgures show the results of the numerical program. Figure El is the 
BPS case where /9 = 1, Figure Elis the type I case with f3 = 1/16 and Figure Elis the 
type II case with jS = 16. In every case we present the prohles for small n = 100 and 
large n. At n = 25, 000 the step function for x(r) is evident. Note that in the large 
n limit the information of (3 is almost lost and the prohles are all the same. The case 
n = 100 is interesting because we can see the difference among BPS, type I and type 
II while approaching the wall vortex limit. 

As another check we now evaluate the tension of the vortex. If we directly insert 
the ansatz (1221) into the energy density derived from the Lagrangian dnD, we obtain; 


Tv{n) = 271 rdr 


1 

L2e2 V 


/nA' 


n 


+ X(1 - + P + r (q) 


(3.5) 


Now we want to verify the conjecture in its hrst formulation ( 1221 ), that is by increasing 
n without changing the parameters of the theory. We introduce a new quantity T (n) 
dehned to be the ratio between the real tension and the BPS tension: 

Tv{n) 


T{n) = 


(3.6) 


7bps (’^) 

For this new quantity the conjecture (I22D is simply lim^^oo 7*(n) = 1. 

Figure El shows the three plots of Tin) for (3 = 1/16,1,16, respectively. The 
quantity T{n) has also an important physical meaning: it is the tension per unit of 
flux carried by the vortex. This means that from the sign of the derivative dT / dn, 
we can read if there is attraction or repulsion between vortices. Figure El is consistent 
with the ordinary expectation of type I and type II superconducting vortices. If we 
take for example (3 = 1/16 (type I vortices) the derivative dT/dn is negative and 
this means that there is attraction. For /? = 16 the function T grows up to 1 and 
this means that there is repulsion (type II vortices). Our theory predicts also the 
way T{n) approaches 1 at inhnity. The radius of the vortex in the large n limit is 
i?MiT y/n while the thickness of the wall A^/ remains constant. This mean that 
the deviation from the “perfect” wall vortex is of order 1/y/n. In Figure El we show 
also the hts T (n) ~ 1 + const /^/n and they perfectly agree with the numerical data. 
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n=8.5x10^, p=1/16 


Figure 2: x(r) and A{r) profiles for (3 = 1/16 and n = 100, n = 85, 000. 



Figure 3: x(r) and A{r) profiles for /3 = 1 and n = 100, n = 25,000. 



Figure 4: x(r) and A(r) profiles for /3 = 16 and n = 100, n = 25, 000. 
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1.8 


Tension normalized to BPS for various values of p 


p=1/16 + 

P=1 
p=16 



0x10° 2x10^ 4x10^ 6x10^ 8x10^ 1x10° 


n 


Figure 5: Plot of the normalized tension T as function of n for the three cases f3 = 1/16,1,16. 
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Figure 6: Plot of the normalized tension T in the log-log graph. 

3.2 Degenerate Vacua and the SLAG bag regime 

Now we consider the case of degenerate vacua when the potential is like that of Figure 
[71 Since the Coulomb vacuum energy density vanishes, there is no MIT regime in the 
large n limit. The conjecture (1231) must then be replaced by 

lim Tv{n) = TslacW • (3.7) 

n—>cxD 

If we want to test the conjecture with our numerical code, it is necessary to hnd 
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Figure 7: The Coulomb phase and the Higgs phase are the two degenerate vacua of the potential. 


a rescaling of the parameters so that the radius of the vortex remains hxed while the 
thickness of the wall converges to zero. Since the large n limit is a SLAG bag regime, 
we must properly modify all the analysis that we have done in subsection 12.21 This 
has already been done in P] so we present only a brief description. The equation 
(ITT^ must now be substituted with 


V(,) ^ (1) . 


(3.8) 


The differential equations (12.1 :-iH2.1411 are the same and the three relevant parameters 
are 

n , a = , b = 2e^qo^ . (3.9) 


2go2 

The radius of the SLAG bag (ESD can be rewritten in terms of the three relevant 
parameters 


-Rslac 


n 


2/3 


(3.10) 


ai/6),i/3 ■ 

We can now give the new formulation of the conjecture. 

Consider the succession of parameters and bn = n^^^bi and call the 

solution of \2.1di) and \2.1f\l with the vortex boundary conditions, and 
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An,a„,bn{i^) ■ In the limit n ^ oo 


lim Xn,ar,,bAr) ^Onir- R^) 
hm An,a„,bAr) ^ ^ 


0 < r < Rv , 
r > Rv . 


(3,11) 


This limit has been chosen so that the radius of the vortex remains constant Ry ~ 

and also the ratio an/bn remains constant. The ratio 
a/h disappears in the wall limit and is only related to the shape of the limiting 
functions Xn,a„,b„(?^) and 

For the computation we take the simplest potential with the Coulomb and Higgs 
degenerate vacua: 

V(x)=x^(x^-1)^ (3.12) 

In Figure IHl we have plotted the tension per unit of flux T{n)/n for three different 
values of a/b. As predicted by the theory all the lines converge to T/n oc as n 

becomes large. 


Tension over n of degenerate vacua with various values of a/b 



Figure 8: Tension per unit of flux for three different values of a/b. At large n, T cc n as 
predicted by the theory. 


In Figure ini we have instead plotted T{n)/n‘^R. In this way is more easy to see 
the SLAG regime since at infinity all the three lines approaches a constant. In this 
plot we can also see the first order deviation from the “perfect” SLAG bag. Since the 
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Tension over n‘^'^ of degenerate vacua with various values of a/b 



OxlO"" 1x10*^ 2x10 3x10*" 4x10 SxlO*" 6x10“ 7x10“ 8x10 9x10“ 


Figure 9: Tension divided by for three different values of a/b. At large n they approach a 
constant with an error of order l/n?'^^. 

radius -Rslac ot while the thickness of the wall remains constant, we expect a 
deviation of order The plots with corresponding hts conhrm this expectation. 

An interesting question is if there exists a kind of BPS notion also for the degen¬ 
erate vacua potential. Maybe there do exist some potential and some particular value 
of a/6 for which the tension not only approaches Tslac as n goes to inhnity but is 
exactly equal to it for all values of n. 

3.3 SLAC/MIT phase transition 

Another thing to be checked is the presence of a window, namely the SLAG regime, 
when the Coulomb vacuum is metastable and the peak vq is much greater than the 
energy density ^o- To obtain this regime we need a potential like the second of Figure 
Qwith Vq ^ £o- The theoretical analysis predicts that in the windows of parameters 
(1213) the tension and the radius scale like eq. (j2ini). 

The simplest way to obtain such a potential is to add an opportune y® interaction: 

V(x) = (ax2 + l)(x^-l)^ (3.13) 

The conditions V(0) = 1 and V(l) = 0 leave one free parameter a that essentially 
measures the height of the peak of the potential. When a is great enough the peak 
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is at X = 1/v^ and has the value Vq = Aa/27. 

We need quite a big a to detect the window (EUD- In Figure Uni we have plotted the 
profile x(^) for o’ = 500. For this computation we have used the rescaling = noi, 
bn = nbi- In this rescaling the large n limit should be a step function with constant 
radius i?MiT = 1- What we should see, to confirm the theory, is that the x(r) prohle 
becomes a step function much before reaching the MIT regime. The radius of the 
wall vortex should be governed by the SLAG law and should increase as When 
the radius reaches 1 there should be a phase transition between the SLAG regime and 
the MIT regime where the radius remains constant as 1. 

Using the rescaling a„ = nai, bn = nbi the numerical computation cannot reach 
high enough n to reach the phase transition. Nevertheless it is possible to detect the 
presence of the SLAG regime just by looking at Figure^] The profile x(^) becomes 
a step function at a radius much lower than 1 and then the radius increases towards 

1 . 

Profile functions for x with p=1 and o=500 Profile functions for x with p=1/16 and cr=500 



Figure 10: Profile yfr) for the potential l.'t.l ;tll with cr = 500. The two graphs correspond to /3 = 1 
and /3 = 1/16, respectively. 

To see the SLAG/MIT phase transition we have to use another strategy. We use 
the rescaling adapted to the degenerate vacua situation bn = n^^^bi and 

we take the potential to be 

v(x) = (d + l)(y - 1)" (3.14) 

cr 

SO that a high value of a corresponds to a small perturbation of the degenerate vacua 
potential dsm. With this strategy it will be possible to detect the phase transition. 
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In Figure [TT] we have plotted the tension for the degenerate vacua potential ()d.l2|l 
and the perturbed degenerate vacua potential (IXTl) with (j = 100. The tensions are 
the same up to n ~ 10^ and then the perturbed potential starts to deviate from the 
unperturbed one. Unfortunately the computation stops here and we cannot see that 
it enters in the MIT bag phase. Anyway the phase transition is clearly visible from 
the Figure. 


0.1 


Phase transition between the SLAG and MIT regime 

perturbed with 0 = 1 00 + 

unperturbed 







Figure 11; In this graph it is possible to see the phase transition between the SLAG regime and 
the MIT regime. The green points correspond to the degenerate vacua potential that is always in 
the SLAG regime. The red points correspond to the perturbed potential Id. 1411 with ct = 100. At 
n ~ 10^ there is a second order phase transition. 


4 Conclusion and Further Developments 


One of the most surprising aspects of the large n limit (EH), is that the tension 
depends only on the value of the potential at zero U(0) = eo- If we have two different 
potentials with the same eo, the large n limit of the tension is the same. 


Now we are going to argue something more. If we take two different potentials 
Vi(g) and 14 (g) with the same Eq and the same go, not only the tensions are the same 
in the large n limit, but also are the prohle functions. 


The energy density derived from the Lagrangian m is the following 


Ti[A, q] = 27i f rdr ^ 


(4.1) 
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where the index i = 1,2 refers to the potentials Vi and V 2 , respectively. The energy 
density must be regarded as a functional: it takes the two profile function A{r) and 
q{r) and it gives out a number. The minimization of Ti[A,q] gives the prohles and 
the tension of the vortex. Now call Ai[r) and qi{r) the profiles obtained by the 
minimization procedure for the potential Vi, and put them into the functional of the 
second potential. What we obtain is 


T2[Ai, qi] — Ti[Ai, qi] + 27r 


rdr 


V2{qi 



(4.2) 


where we have simply added and subtracted Vi. In the large n limit the profile qi{r) 
becomes a step function and, since the two potentials have the same eo and go, the 
extra piece 27r J rdr[ 1 / 2 ( 51 ) — l/i(gi)] vanishes. This means that the functions Ai{r) 
and 51 (r), being the minima of the functional Ti, become also approximately the 
minima of the functional T 2 . 

If the potential 1 / 1 ( 5 ) ^ 2 ( 5 ) have the same Eq and the same 50 , the profile 

functions converge in the large n limit 


lim II 51 - 52 II = 0 , lim ||Ai - A 2 II = 0 , (4.3) 

n—^■OO 72—^■oo 

where the distance between functions is measured with the If metric. If the conjecture 



V{\q\) 


A 

£0 



Figure 12: In the first sketch we have two different potentials with the same Eq and the same qq. 
In the second sketch we have two potentials with the same Eq but different qo- 

is true we have also the following result. If the ratio between the zero energy density 
Eq and the vev 50 is the same of the BPS potential, that is Eq = eqo/2, whatever the 
shape of the potential, we recover supersymmetry in the large n limit. It is in fact 
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known that the Abelian-Higgs model with the BPS potential arise from the bosonic 
Lagrangian of SQED with zero mass and a Fayet-Ilyopoulos term In this case 
the “miracle” of the proportionality between the tension and the charge finds its 
explanation in the central charge of the supersymmetry algebra |E|.3 

In the present paper we have hnaly found a convincing proof that multi-vortices 
become bags in the large n limit. This fact is very general and applies also in the 
multi-monopole case ^Ij. It will be interesting to see if the bag mechanism works 
also in more general theories that contain solitons: nonabelian vortices PI. semi-local 
strings PI, Chern-Simon theories PI and noncommutative field theories m- 
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